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ABSTRACT 

We explore the dynamical behaviour of cosmological models involving a scalar field (with an exponen¬ 
tial potential and a canonical kinetic term) and a matter fluid with spatial curvature included in the 
equations of motion. Using appropriately defined parameters to describe the evolution of the scalar 
field energy in this situation, we elucidate the character of two fixed points that are not present in the 
case without curvature. We also analyse the evolution of the effective equation-of-state parameter for 
different initial values of the curvature. 


1. INTRODUCTION 

The discovery, based on the observed behaviour of Type la Supernovae [1, 2], that the expansion of the Universe 
appears to be accelerating, generated enormous theoretical interest in finding a suitable framework to account for this 
phenomenon. Other independent observational data such as the Cosmic Microwave Background radiation (CMB) [3, 4] 
and Baryonic Acoustic Oscillations (BAO) [5, 6] have subsequently confirmed that reconciling a standard Friedmann- 
Lemaitre-Robertson-Walker model based on Einstein’s General Theory of Relativity requires roughly 70% of the 
average energy density today to be in the form of an exotic fluid whose equation-of-state parameter w = P/p is close 
to —1. 

Although over fifteen years have passed since the original discovery, we still lack a compelling theoretical model to 
explain this so-called dark energy. The three major theories that attempt to do this all rely in one way or another on 
modification of the Einstein Equations, which we use in the form 

Rpu - Rg^v/2 = 8ttGT^, (1) 

in which properties of space-time appear on the left-hand side and the matter-energy contents of the Universe on the 
right. 

The obvious phenomenological possibility is to add a constant term A g^ v to Einstein’s Equations. This seems 
appealing because of its simplicity and that when included in the energy-momentum tensor on the right-hand-side of 
the Einstein Equations, it leads to a clear prediction of the effective equation of state, namely that w = — 1 exactly. 
There are, however, two major concerns with this picture. The first is that the magnitude of the vacuum energy 
associated with a A-terrn is out of line with the value expected from summation of zero-point energies of quantum 
fields up to some cut-off scale. This discrepancy is 120 orders of magnitude if one chooses a cut-off at the Planck 
energy, but even if some mechanism imposes a cut-off at for example a QCD energy scale, the problem is alleviated 
but not entirely solved [7]. 

The second problem is that there seems no natural explanation of why the energy density associated with the vacuum 
energy should be within a factor of a few of the present matter density or, in other words, why the expansion of the 
Universe should have began to accelerate so very recently in cosmic history. These are called the fine-tuning problem 
and the coincidence problem, respectively [8]. Although such statements clearly depend on some choice of measure, 
related in this case to the probability distribution of the value of A, they must also take account of anthropic selection 
effects [9]. 

A second class of models that might explain dark energy are those based on some form of modification of Einstein’s 
theory of gravity. The usual approach in such models is that instead of starting from the standard Einstein-Hilbert 
action (which leads to the Einstein Equations), one considers additional terms in the action [10]. The most straight¬ 
forward modification of General Relativity is to replace the Ricci scalar, R : in the Einstein-Hilbert action, by ther 
function of this scalar, usually called f{R). An example of this model is f{R) = R + aR 2 , which was in fact one 
of the first models of inflation, proposed by Starobinsky. The idea of using f(R) for late-time acceleration was first 
suggested in [11] and examples of viable models of this type are proposed in [12, 13, 14]. In a more general case, 
one can include an arbitrary function of R, R^R^ and Rp,upaR PlVpa in the action. The Gauss-Bonnet combination 
R 2 — 4R^R^ + Rp,v pa R p ' Jp ' J [15] is particularly widely explored in the literature. Another, more complicated, class 
of modified gravity models include scalar-tensor theories, where Ricci scalar R and a scalar field (j) are coupled (an 
example of this is Brans-Dicke theory [16]), and DGP (Dvali, Gabadadze and Porrati) braneworld model [17]. In 
the latter, particles are confined to a 3-dimensional brane, embedded in a 5 (or more)-dimensional space-time with an 
infinite extra dimension. Standard 4D gravity is recovered at small distances, but at larger scales gravity is weakened, 
because its energy is essentially getting lost to the additional dimension. 

A third class of possible models include modifications to the standard form of matter on the right-hand-side of the 
Einstein Equations, designed to generate a negative effective pressure. Among such models are quintessence, fc-essence, 
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coupled dark energy and the generalised Chaplygin gas (for reviews see [18, 19, 8, 20]). Quintessence, which will be 
studied in much greater detail in the rest of this paper, represents the idea that accelerated expansion is driven by a 
canonical scalar field (f> [21, 22, 23]. The most important consequence of this is that we now have a dynamical equation 
of state, rather than a constant [24, 20]. Please note that this is also the case in some of the modified-gravity theories 
discussed in the paragraph above. Quintessence models are usually divided into two types on the basis of the form 
of potential that drives the scalar field dynamics. The first class contains models in which w gradually decreases to 
w = — 1; these models are called “freezing models” [25]. The forms of potential needed for this kind of behaviour were 
studied in [26, 27]. In this scenario, the field energy density does not necessarily need to be negligible at the radiation 
epoch (this is the case with the cosmological constant). The other option is the class of “thawing models” [25] in which 
the field equation-of-state parameter is close to w = — 1 initially, but at late times deviates from this value. 

Theories that involve non-canonical kinetic terms in the Lagrangian are called fc-essence [28, 29, 30]. The idea is 
that inclusion of non-canonical terms results in cosmic acceleration even without the field potential, as was shown in 
[31]. A number of different scenarios were proposed, to name a few: Low energy effective string theory [32], Ghost 
condensate [33], Tachyon field [34], Dirac-Born-Infeld (DBI) theory [35, 36]. 

Furthermore, there are suggestions that dark energy and dark matter are coupled [37, 38, 39], that a single fluid (e.g. 
generalised Chaplygin gas) is responsible for dark energy and dark matter [40, 41], and finally unified quintessence 
and inflation theories, namely quintessential inflation [42, 43]. In [44] the authors consider models of quintessence 
interacting with dark or baryonic matter. 

There are also attempts to explain the apparent cosmic acceleration by means of inhomogeneities in the matter 
distribution in other words by breaking the assumption that cosmological space-time is described by a FLRW metric; 
for a review see [45]. One specific example is the suggestion that we live in the middle of a huge underdensity (a 
“void”) and we interpret the expansion of its surroundings as an overall cosmic acceleration [46, 47], although this 
seems to be in conflict with observations, see for example [48]. Another approach relies upon the back-reaction of 
cosmological perturbations [49, 50], which may be able to explain acceleration without dark energy; some work along 
these lines related to this paper can be found in [51, 52]. 

Given this plethora of possible models it is important to undertake a rigorous systematic study of their dynamical 
properties, in order to understand and classify the wide range of behaviours they may exhibit. In particular, focusing 
on the fixed points of their evolution will allow us to tackle the difficult question of what can be considered to be 
“natural” behaviour in a given scenario. Since the fixed points attract trajectories from a wider parameter space 
around them, the configuration of the system can be determined by its dynamical evolution in the fixed point rather 
than in the initial conditions. We concentrate on a specific class of possibilities by studying the dynamical evolution of 
quintessence-type models, extending previous work (described in detail below) by including spatial curvature, which 
is usually neglected in such analyses. Although the observational evidence at the moment points to a universe which 
is (nearly) spatially flat this conclusion is based on a model with restricted set of parameters. In the framework of a 
more general model the constraints on CLk can be weaker [53]. It therefore remains important to establish whether 
the inclusion of curvature leads to any qualitative changes in the dynamics. This issue has been addressed before [54], 
but in this paper we use a different parametrisation of the equations of motion, which makes the connection between 
the curvature and dynamics more explicit. The aim of our analysis is to ask the question whether it is possible, via 
this relatively simple generalisation of the quintessence scenario, to generate an attractor solution that corresponds to 
a value of w ^ —1, which would distinguish it from many of the other models listed above. 

The paper is organized as follows. In Sec. 2 we explore the dynamics of a single scalar field with an exponential 
potential and a canonical kinetic term under the flat FLRW metric. This is a well-studied system [55, 56] but one 
which repays further analysis. We find the fixed points of this system and investigate values of important dynamical 
parameters at these fixed points; this section follows the analysis that was done in [24]. The further details of this 
analysis are given in the Appendix. In Sec. 3.1 we generalize this approach by introducing another variable that 
corresponds to spatial curvature and analyze how the dynamics are affected as a consequence. We find two new fixed 
points and explore evolution of equation-of-state parameter for different initial values of curvature. For a parallel 
discussion see [57]. 


2. DYNAMICAL ANALYSIS 
2.1. Background 

The action of a canonical scalar field (j> is given by 

S = j d^XsJ—g 


\m^R - - V(0) + Cm 


( 2 ) 


where Mp\ represents the reduced Plank mass: Mp\ = l/y/8nG = \/yfk. The matter Lagrangian Cm has two 
contributions: non-relativistic matter with equation of state w m = Pm/Pm = 0 and radiation with equation of state 
w T = 1/3. Varying this action with respect to metric and applying the action principle gives the Einstein Equations 
with the energy momentum tensor: 


= d^(j)d v (t> - g 


^d p (\)d p (/>- V{4>) 


( 3 ) 
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The assumption of spatial homogeneity and isotropy allows us to adopt the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metric: 

dr + r 2 dd 2 + r 2 sin e 2 d(t> 2 


ds 2 = —dt~ + a(tY 


1 - Kr 2 


( 4 ) 


where a(t) is time-dependent scale factor and K is a constant describing the spatial curvature. The universe is open 
if K < 0, flat if K = 0 and closed if K > 0. With the use of this metric, the Einstein Equations become the modified 
Friedmann Equations: 


3 Mp\H 2 — -(f ) 2 + V (</>) + Pm + Pr — 3Mpj — 

2 cr 

2iVfpjij = — </> 2 — (1 + w m )p m — (1 + w r )p r + 2-Mpj — 


(5a) 


(5b) 


where the Hubble parameter is defined in the usual manner as H = a/a. The requirement of energy-momentum 
conservation demands = 0. For a homogeneous and isotropic universe the energy-momentum tensor becomes 

symmetric, i.e. T= diag(— p, P, P , P ) and we thus obtain the continuity equation p + 3 H(p + P) = 0. This equation 
is obeyed separately by the matter, radiation and scalar field as long as there is no coupling between these components. 
Applying the same condition to the energy-momentum tensor of the scalar field (3) we get 

4> + 3 H/j) + V t4> = 0 , (6) 

where V^, is the derivative of the potential with respect to the field 4>. Comparing this to the general continuity equation 
we see that the effective energy density and pressure of the scalar field are p<f> = </> 2 /2 + V(<j>) and P$ = <^ 2 /2 — V(<j>), 
respectively. The equation of state for the scalar field is then given by 

<fi 2 /2 - V(4>) 


p< 

Wtjy = - = 


P<t> 0 2 /2 + V(<t>) 


( 7 ) 


Equation (6) is a dynamical equation for the evolution of the scalar field; in order to solve for the dynamics of the 
system, one must solve this equation simultaneously with the Friedmann Equations (5) and the continuity equation. 


2.2. Two dynamical variables 

It is convenient to define a new set of dimensionless variables: 


x = - , y = 

V6M Pl H 

because the energy density of the scalar field can be expressed as: 

n P 0 2,2 

^ ~ 3Mp[i? 2 " + V 


VW) 

V3M Pl H 


and the equation of state for the scalar field reads: 


2 2 
2 T — y z 

W6 = -9 

x z + y 


2 ' 


Additionally defining combined matter and radiation energy density parameter as: 

Pm , Pr 


— 


3M 2 X H 2 


3M 2 X H 2 


( 8 ) 


(9) 


( 10 ) 


(ID 


and assuming that spatial curvature is zero, the first Friedmann Equation (5a) simplifies as: 

1 = x 2 + y 2 + (12) 

To study the evolution of the scalar field we can take derivatives of x and y with respect to the number of e-foldings 
N = In a, anticipating that such a system can display accelerated exponential expansion. For an exponential potential 
V = H oe -W/M pi we get . 


= \ZI Ay2 " \ x ( Wm + v 2 - 1) ^ 2:2 + v 2 +1) 

= -\j\ Xx V - \y ( W M (x 2 + y 2 -l) -x 2 + y 2 . 


(13a) 


(13b) 


In order to analyse the dynamical system (13) we first look at the fixed/critical points for different values of the 
parameters A (not to be confused with the cosmological constant A) and (defined as wm = (Pm®m+Pr® r )/(Pm+Pr))- 
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From the constraint (12), it follows that x 2 + y 2 < 1 because Om is always positive for the form of potential we use. 
Every solution with non-zero y can be positive or negative in y , since y 2 ~ V. The part of parameter space in which 
y is negative corresponds to a contracting universe. Because the phase plane is symmetric with respect to the x-axis, 
we will only consider its upper part. This implies that trajectories in phase plane are limited to the upper half of the 
unit disc. Fixed points are defined as: 




(14) 


There are five such critical points; their coordinates and values of the associated physical parameters are listed in 
Table (1). 


TABLE 1 

Coordinates and properties of fixed points. 



X * 

y* 

'R’eff 



A 

0 

0 

WM 

0 

1 

Bi 

1 

0 

1 

1 

0 

P2 

-1 

0 

1 

1 

0 

D 

A/%/6 

%/l - A 2 /6 

A 2 /3 — 1 

1 

0 

E 

V3/2(w m + 1)/A 

•\/ 3/2(1 — 


3(ip m + 1)/A 2 

1 - 3 (w m + 1)/A 2 


The value of the parameter A can consequently be divided into three qualitatively different cases: 


• (A < a/3(u;m + 1)) All trajectories are drawn to the point D which is a stable attractor. 

• (a/3(wm + 1) < A < a/6) Point E becomes a spiral attractor. D is still present. 

• (a/6 < A ) D is no longer defined and Bi becomes a saddle point. 

The solutions of Equations (13) for these three cases are illustrated in figure (1). We have defined ®m to be a 
combination of pressureless dust (with k;m = 0) and radiation (him = 1/3) so its values can only lie between these two 
values but in the most general case (including more exotic equations of state) its value could lie outside this range. 
The specific cases shown in the figures (1) all have kjm = 0. The value of this parameter only slightly affects the 
position of the fixed point E (and none others). For higher values of «jm points E and D merge at higher A. For a 
limiting case ium = 1 (that extends beyond our analysis) point E merges with D exactly at the point Bi so it is never 
inside the half disc. 



Figure 1. Phase planes for three qualitatively different cases of A and wm = 0 as discussed in the text. 
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3. THE ROLE OF CURVATURE 
3.1. Three dynamical variables 


In this subsection we introduce spatial curvature into dynamical system, while still keeping matter and radiation as 
one variable and considering a scalar field described as an exponential potential. We define new variables to construct 
a generalization of the two-dimensional phase plane discussed in the previous subsection into a three-dimensional phase 
space: 


0 

y/6M Pl H ’ 


? JvW) 

V y/3MpiH ’ 


K 

a 2 H 2 ' 


(15) 


Definitions of f2^, w$, Dm and % are identical to those in the previous subsection. The first Friedmann Equation is 
now expressed as: 

1 = x 2 + y 2 - z + D m . (16) 


Note that z is not squared in this definition. If 2 were defined in such a manner that it was squared in the above 
equation then all solutions with negative K would correspond to solutions with imaginary z. The minus sign in front 
2 has no physical significance; our definition corresponds to Qk = —z. With the introduction of the curvature, ratio 
of the Friedmann Equations becomes: 


\ ( w m (x 2 + y 2 - z-l) - x 2 +y 2 -l) - *z 

and the effective equation-of-state parameter can be expressed as: 

P(j> + Pm x 2 - y 2 + w M (l - x 2 - y 2 + z) 

Weff = - = -y—- 

p(j) + Pm 1 + z 

The evolution of the dynamical system is now described by a system of three coupled differential equations: 


(17) 


(18) 


^ = \J\ X V 2 ~ \ x ( 3w m (x 2 + y 2 - z - 1) - 3x 2 + 3y 2 ~ z + 3) (19a) 

^y (-3ic m {x 2 + y 2 ~ z — l) + 3x 2 - V6\x - 3 y 2 + z + 3^ (19b) 

cLz 

—T = 2 (— 3 k;m (x 2 + y 2 - z - l) + 3a: 2 - 3 y 2 + 2 : + l) . (19c) 

dN ' ' 

3.2. Phase space analysis, fixed points and their stability 

In this case under consideration the constraint x 2 + y 2 < 1 holds only in the plane z = 0. The general constraint is 
x 2 + y 2 — z < 1 and it defines a parabolic surface in the space of parameters x, y, z. All trajectories must lie above 
this parabolic surface; they are confined either to the plane z = 0, where the curvature is zero, or to the part of the 
space where z > 0 and K has a positive sign, or to the finite part of the space between the parabolic surface and the 
plane z = 0 with a negative ^-component and thus a negative curvature parameter K. Because this parameter can not 
change sign in our model, no trajectories can cross from 2 > 0 to 2 < 0 or vice-versa. Figures (2) show solutions for 
four different case of parameter A. Again, we ignore u>m by setting it to zero. 

Inspection of the fixed points in this new system shows that all fixed points we had before are preserved, when the 
newly-added additional 2 is fixed at z = 0. In addition, however, there are now two new fixed points: C and F. The 
eigenvalues corresponding to these points are obtained by extending the matrix (8) to describe three variables; the 
additional dimension of the parameter space requires that there will be an additional eigenvalue. If that eigenvalue is 
positive, then the fixed point is repulsive in the new direction and if it is negative, the fixed point is attractive. For a 
point to be an attractor, all the eigenvalues must be negative at that point. If at least one is positive, then trajectories 
will be drawn away in at least one direction. 


Point A: is the matter dominated solution where all energy density is in Dm- At this point the only contribution to 
the effective equation of state parameter comes from the matter sector, so w e g = Wm- The eigenvalues of this 
point are [3(«)m — 1)/2,3(ium + l)/2, (1 + 3u>m)] so the point is a saddle for all values of A, as long as % < 1- 
It is attractive for the trajectories along the a;-axis and repulsive for others. In the limit case of initial potential 
being exactly 0 (the trajectories that start exactly on the a;-axis) this point is an attractor. 

Points B x and B 2 : represent solutions in which the universe is dominated by the kinetic energy of the scalar field: 

= 1. Here the effective equation-of-state parameter is constant and the scale factor behaves as a oc t 2 / 3 . 
Eigenvalues for stability are [3(1 — %), i/3/2A =F 3,4] with the minus sign for B x and plus for B 2 . This means 
that B 2 is always repulsive and Bi repulsive for A < \ff> and a saddle for A > \/6. 
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Point C: is new compared to the two-dimensional case and corresponds to a curvature-dominated solution; all the 
energy density resides in f7jc- This solution is trivial, because it does not involve any matter or the scalar field. 
The point is a saddle for all values of parameters; its eigenvalues are [—2,1, —(1 + 3 wm)]- 

Point D: is a solution where the universe is dominated by the scalar held. All energy density is in 17^ = x 2 + y 2 = 1, 
with Om = 0. The effective parameter of state is w e s = A 2 /3 — 1, so the universe is accelerating for A < y/2. 
In the limit case A —> 0, this solution corresponds to the de Sitter expansion dominated by the cosmological 
constant. This fixed point lies in the z = 0 plane, at the edge of the half disc and it moves from (a: = 0,2/ = 1) 
for A = 0 to (x = l,y = 0) for A = y/6. After this value of A it is not defined anymore. The eigenvalues are 
[(A 2 — 6)/2, A 2 — 3 (u>m — 1), A 2 — 2] which means that the point is attractive for all trajectories with z < 0 and 
some with z > 0 for A < y/2. It is also an attractor for trajectories in the plane z = 0 for \[2 < A < \/3(wm + 1), 
and a saddle for + 1) < A < y/6. The case of A < a/ 2 is especially interesting because in this setup even 

some of the trajectories that start with a positive amount of curvature end up in the inflationary solution (most 
of the trajectories which start with the positive curvature evolve to a state with infinite K, i.e. the universe 
collapses). This demonstrates that “closed” does not imply “finite”. 

Point E: is the so-called “tracking” solution. Both 17^ and 17 m are between 0 and 1, and the effective equation-of-state 
parameter matches that of the matter: w e s = Wm- This means that the universe expands as if it were matter 
dominated and there is no accelerated expansion. This point is, however, interesting because the scalar field is 
still present; some of the energy density in stored in f7^. This fixed point lies in the z = 0 plane where both its 
coordinates are infinite for A = 0. Since that is outside of the half disc, it is not relevant for this analysis. The 
point only becomes relevant as A becomes larger than y / 3(wM + 1)- At this value E crosses D and enters the 
half disc. Analysis of the stability shows that for A > y/3(iCM + 1) two of the eigenvalues are complex conjugates 
of each other. This means that for the trajectories in the z = 0 plane this point is a stable spiral. The additional 
eigenvalue is always positive so all other trajectories are repelled. 

Point F: The effective equation-of-state parameter in point F is always w e j = —1/3. This corresponds to the scale 
factor being linearly proportional to t. The matter-energy density vanishes here so the point always lies exactly 
on the parabolic surface x 2 + y 2 — z = 1. The corresponding eigenvalues are [—3 wm — 1, — v 7 8A 4 — 3A 6 /A 3 — 
1, a/ 8A 4 — 3A 6 /A 3 — 1] which results in this point being either a saddle, attractor or a stable spiral, depending on 
the value of A. This is the analogue of the tracking solution (point E) for the case where the universe behaves 
as if curvature were the only component. This implies a oc t. 


TABLE 2 

List of fixed points and theirr properties for the three dimensional case. 



X * 

y* 

Z* 

^eff 



Mr 

A 

0 

0 

0 

1PM 

0 

i 

0 

Bi 

i 

0 

0 

i 

1 

0 

0 


-i 

0 

0 

i 

1 

0 

0 

c 

0 

0 

-1 

undefined 

0 

0 

1 

D 

A/V6 

a/1- A 2 /6 

0 

A 2 /3 — 1 

1 

0 

0 

E 

y/3/2(w M + l) 

X 

V 3 /2(l-0 

A 

0 


3(ui m + 1) 

A2 

-I 3(ui m + 1) 

1 P 

0 

F 

a/2/3/A 

2/(v / 3A) 

2/A 2 - 1 

-1/3 

2/A 2 

0 

1 - 2/A 2 
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A=1.2 


A=1.6 



Figure 2. Phase spaces for four different A and wm = 0. Valid phase space is limited by the parabolical surface: all valid trajectories 
lie above it. Trajectories in the x — y plane (where z = 0) are shown in blue. Green trajectories lie between this plane and the parabolic 
surface. Since z = —FIk — K/(a?H 2 ), green trajectories represent part of the phase space with negative curvature. Red trajectories are 
bounded to the part of the space where z > 0 which corresponds to K > 0. Point E lies under the parabolic surface for the first two cases, 
so it’s not shown in the plots. 






















Four qualitatively different cases of A are distinguishable: 

• (A < y/2) There are 6 fixed points (all except E, which is outside of allowed region). Point D is the attractor for 
all of 2 < 0 part of the space and some trajectories in z > 0. F has positive ^-component and is a saddle. 

• (\/2 < A < y / 3(wM + 1)) There are still the same 6 fixed points but the 2 -component of F is now negative. This 
point is a stable attractor for all trajectories with z < 0. For trajectories that lie in the plane 2 = 0 the attractor 
is point D, and trajectories with 2 > 0 don’t converge. 

• (y/Z(wM + 1) < A < \/6) Point D becomes a saddle, E enters the allowed region and becomes attractive for 
trajectories in the plane, while trajectories with z > 0 diverge. Two of the eigenvalues of point F are now 
complex conjugates of each other while the third is negative, so this point is a stable spiral in three dimensions. 

• (y/6 < A) At A = y/6 The significance of this last case is that point D merges with Bi and disappears. Bi is 
now a saddle point in the plane. 


3.3. Evolution of the Equation of State 

For the same cases of A we plot trajectories that correspond to initial conditions that give flat, open and closed universe 
models. These are shown in figures (3). Note that we used the same initial conditions for all four cases. In the first 
case all curves converge. For all other values of A the equation-of-state parameter diverges for I\ > 0. For K < 0 
trajectories in phase space are drawn towards point F, where w e g = —1/3. For K = 0 they can either end up in the 
point D (first two values of A) or E (third and fourth case of A). 



g 

£ 


— K >0 

— K =0 

— K <0 


0 5 


10 15 20 

N 





Figure 3. Evolution of equation-of-state parameter over 20 e-foldings (N) for four different cases of A. 
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4. CONCLUSIONS 

We explored a single-field quintessence model with exponential potential and canonical kinetic term in presence of a 
FLRW metric and overall (positive or negative) spatial curvature. In the Friedmann Equations, the term with spatial 
curvature is dynamical, so we introduced a new variable and generalized two-dimensional case that holds for a flat 
geometry to three dimensions. 

In comparison to the model without curvature there are two new fixed points. These have been noticed before (e.g. 
[56]) but their behaviour is clearer when viewed in terms of the parametrization we use. One of them is trivial, in 
that it corresponds to a universe where entire energy density is dominated by the curvature. This fixed point is not 
attractive. The other fixed point corresponds to a universe where energy density is a combination of curvature and 
scalar field (but no matter). The ratio between them depends on the parameter of the exponential potential. This 
fixed point is interesting, because it is attractive for all trajectories with negative curvature term for A > \/2. We have 
thus established that there is a natural configuration of this system that corresponds to a specific value of w. Although 
we know that the energy density of the Universe today is almost critical (so there is very little or no curvature), this 
fixed point might still be interesting for the cases where A is just above y/2 and so energy density is dominated by the 
scalar field. 

The only problem with relating this solution to our observed reality is that equation-of-state parameter in this point 
is —1/3, but we know that the present value is very close to —1. Nonetheless, the attractor solution D is an accelerating 
solution for flat, negatively curved and some positively curved cases. The corresponding equation-of-state parameter 
is close to the value w = — 1 in the limit A —► 0, but for any value of A < x/2 it corresponds to w < — 1/3 so this is at 
least qualitatively applicable to dark energy. 

We have only investigated one particular and rather simple model, but this sort of curvature inclusive analysis is 
applicable to other, potentially more complicated models that might exhibit a much richer dynamical interplay. Cases 
which might prove amenable to further study on these lines would be scalar fields with a more general potential [58], 
non-canonical terms, with multiple scalar fields, with globally anisotropic metrics (i.e. the Bianchi models) [59], e.g. 
the Kantowski-Sachs model [60], and models based on exact inhomogeneous cosmologies. Such models would of course 
introduce more than one extra parameter, so the resulting phase portraits would involve more than three dimensions 
and their analysis would entail considerably greater complexity. 
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APPENDIX 


Two dynamical variables and stability of fixed points 

After introducing new variables as in Equation (8) for a general potential V((/>) and using the Klein-Gordon Equation 
(6), the derivatives with respect to N are: 


dx = <j> / 3 H_\ _ V/ft 

dN y/6MpiH y H 2 J s/6Mp\H 2 

dy _ <j> _ W H 

dN ~ 2v / 3MpiV / Wff 2 V3M Pl H H 2 ' 

Defining parameter A as A = —MpiV^/V is especially convenient because this parameter is constant for the exponential 
potential. Using this potential and eliminating Dm with the use of constraint (12), the ratio of the Friedmann Equations 
H/H 2 can be expressed, in terms of x and y, as: 


(la) 

(lb) 


H 

H 2 


I (w M ( X 2 + y 2 - l) - X 2 + y 2 - l) , 


( 2 ) 


where we define wim to be the combined equation-of-state parameter for a fluid component taken to describe both 
matter and radiation: 


ui M 


PmW m + p r W T 

Pm + Pi 


(3) 


Derivatives of x and y are hence expressed, for the exponential potential as (13). This provides a useful way to express 
the evolution of both and f2^: 


= (w^ - 1) ^-A^3(l + + 3(1 + w^j 


(4a) 


% = 3(Q^ - 1 - ui M ), 


(4b) 
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where the prime stands for differentiation with respect to N. Additionally, by defining T = VV^/V 2 ^ we can obtain 
the derivative of A: 

A' = -A 2 (r-l)^3(l + u;*)fi*. (5) 

For the exponential potential A is constant, so X' does not carry any additional information in this case. In the fixed 
points the effective equation-of-state parameter, defined as 


P<t> + P M 2 2, 2 2\ 

Weft = ---= x - y + u> M (l - x -y ) 

P<j> + Pm 

is constant so we can easily integrate Equation (2) to see how scale factor changes with time: 


( 6 ) 


a OC 


( 7 ) 


Note that this only holds if the system is at one of the fixed points. In other cases Equation (2) has to be integrated 
numerically because its right-hand side is not constant. However, it is still sensible to define an equation-of-state 
parameter w,' e ff even outside these points. For iu e ff = —1/3 the scale factor will have a linear dependance on time; 
values larger or smaller than that will give decelerated and accelerated expansion, respectively. 

We investigate the stability of the fixed points by finding eigenvalues of the matrix, evaluated at the fixed point 

f df* d/s 1 


M = 


dx dy 
dfy ply 

dx °y J x,,y. 


( 8 ) 


where ft = di/dN stands for the terms on the right-hand side of Equations (13). A fixed point is stable (attractive) if 
both eigenvalues are negative, non-stable (repulsive) if they are positive and a saddle (attractive in one direction and 
repulsive in others) if they have the opposite sign. If the two eigenvalues are complex conjugates of each other, then 
the fixed point is a stable spiral. 

In the case of two variables there are five such critical points; their coordinates and values of the associated physical 
parameters are listed in Table (1). Their existence and properties depend on the values of A and i»m- 
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